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Abstract 

Effects of inverse triad corrections and (point) holonomy corrections, occuring in 
loop quantum gravity, are considered on the properties of Reissner-Nordstrom black 
holes. Version of inverse triad corrections with unmodified constraint algebra reveal 
the possibility of occurence of three horizons (over a finite range of mass) and also 
show a mass threshold beyond which the inner horizon disappears. For the version 
with modified constraint algebra, coordinate transformations are no longer a good 
symmetry. The covariance property of spacetime is regained by using a quantum 
notion of mapping from phase space to spacetime. The resulting quantum effects in 
both versions of these corrections can be associated with renormalization of either 
mass, charge or wave function. In neither version is the Newton's constant renor- 
malized. (Point) Holonomy corrections are shown to preclude undeformed version of 
constraint algebra as also a static solution, though time independent solutions exist. 
Possible reason for difficulty in constructing a covariant metric for these corrections is 
highlighted. Furthermore, the deformed algebra with holonomy corrections is shown 
to imply signature change. 

1 Introduction 

Symmetries play a fundamental role in physics. Certain symmetries, like the spacetime 
symmetries, lead to familiar conservation laws for energy and momentum. There are also 
gauge symmetries familiar from particle physics (which actually reflect the redundancy of 
the classical formulation of a physical theory) and are of prime importance, especially in the 
context of quantum theory. In a diffeomorphism covariant theory like general relativity, 
these two notions of symmetry become closely related. This is particularly clear when 
contrasting the usual action formulation of general relativity with the canonical Arnowitt- 
Deser-Misner (ADM) formulation yy. 

While the spacetime covariance property of the theory is obvious in the standard ten- 
sorial formulation, in the ADM formulation spacetime symmetry gives way to gauge sym- 
metry in the phase space and is encoded in the form of constraint algebra of the theory. 
Gauge transformations generated by the constraints then correspond to spacetime dif- 
feomorphisms. It then becomes important to ask as to what happens to the standard 
covariance property of the theory in the quantum domain. Does the symmetry survives 
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in the quantum theory or it becomes anomalous or is it just deformed? In loop quantum 
gravity, where the ADM variables are replaced by Ashtekar variables |2], the issue has 
recently seen renewed interest [31 IH |3] . 

In loop quantum cosmology, where the program of loop quantum gravity has found its 
most successful application [HI [3, [8], the above issue trivializes since the diffeomorphism 
constraint is identically equal to zero. Addition of perturbative inhomogeneities makes a 
full quantum analysis difficult but at an effective level, where inverse triad corrections or 
holonomy corrections are included, it shows that the constraint algebra is deformed P]. 

Like cosmological models, spherically symmetric models have always been a very useful 
test bed for theories of quantum gravity and these have been studied within loop quan- 
tum gravity for different applications (see for instance [TOl [HI US])- These models are 
interesting for several reasons. Black hole geometries in classical general relativity are 
mostly spherically symmetric and thus provide an opportunity to study possible effects of 
quantum gravity on issues related to singularity and horizon properties. Modifications to 
horizon properties can, for instance, have consequences for Hawking radiation, particularly 
during the end stage of black hole evaporation. These models are quite versatile and apart 
from allowing an exploration of black hole spacetimes which mostly have time independent 
geometry, they also allow time dependent gravitational collapse scenarios. 

Apart from these practical considerations there is another very important reason for 
the usefulness of these models. As opposed to cosmological models, the constraint algebra 
is non-trivial here and as pointed out in the beginning, it is the constraint algebra which 
encodes the symmetry properties of the theory. These models therefore provide an oppor- 
tunity to study the effects of quantum gravity on the symmetries of the classical theory 
without entailing all the difficulties that a full theory without any symmetry reduction 
would have. Moreover the inhomogeneity is non-perturbative implying (in certain sense) 
a closer association with full theory. Like cosmological models incorporating perturbative 
inhomogeneities, in several of these models, loop quantum gravity inspired corrections lead 
to a deformed constraint algebra [13], [El [15] . 

Implications of the classical hypersurface deformation algebra from the perspective of 
the canonical theory have been analyzed in pTJ and its implication for the Lagrangian of the 
theory in [18] . Recently in [3l [19] the approach of [18] was adopted for the case of modified 
hypersurface deformation algebra and it highlighted several interesting features including 
generic singularity resolution as a direct consequence of deformation of the constraint 
algebra. 

While the utility of such general analysis is obvious, it is also useful to consider specific 
models since they allow explicit solutions to be obtained which can be compared and con- 
trasted with their classical counterparts and can thus provide useful hints for more generic 
constructions. Such an analysis was performed for spherically symmetric vacuum space- 
times and those incorporating perturbative scalar matter in [15] (see [T6] for a summary of 
some of the results). It was found that with deformed algebra familiar spacetime concepts 
(spacetime metric, for instance) need to be reconsidered as also the constructs like the 
black hole horizon, which are based on such notions. 

In this paper we extend the work in |15j and consider minimally coupled Einstein- 
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Maxwell theory with spherically symmetric background. Reduction to spherical symmetry 
in terms of Ashtekar variables has been discussed in [201 l22] and including Maxwell 
field in [23] (also see ^4J). The aim of the present work is to include loop quantum 
gravity corrections at an effective level and to work out the consequences for the Reissner- 
Nordstrom black hole. The main theme would be to maintain the first class nature of the 
constraint algebra (though it can be deformed). 

In [15] only inverse triad corrections arising from loop quantum gravity were considered. 
Here we also consider the effects of (point) holonomy correction and unlike cosmology, 
the matter field can be non-perturbative. Since the details of the theory for these more 
complicated models are still unclear it is useful to explore various models at an effective 
level incorporating different matter fields. Holonomy corrections have so far not been 
analyzed in much detail for spherically reduced models though they have been discussed 
in some detail in [23]. These corrections can lead to some very dramatic effects leading to 
signature change of the metric as has been pointed out recently in [2S1 E] • 

With quantum corrected constraints and constraint algebra it is not guaranteed that a 
consistent set of equations will result. We will show that at the effective level considered 
here, it is possible to obtain consistent equations for both inverse triad and holonomy 
corrections. Since we want to include the effects of loop quantization which requires a 
Hamiltonian framework, we start in section 2 with the Hamiltonian formulation of the 
classical theory in terms of Ashtekar variables. After performing the constraint analysis 
at the classical level we will move on to consider the effects of inverse triad corrections in 
section 3. As in earlier works we will consider two cases - inverse triad corrections which 
keep the constraint algebra unmodified and those which modify the algebra. For modified 
algebra the familiar spacetime covariance properties do not hold. We recover such a notion 
by providing a quantum version of mapping from phase space to spacetime metric. 

In section 4 we will consider the inclusion of corrections due to point holonomies cor- 
responding to spherically symmetric directions. Here also we will consider two cases - (1) 
phase space independent holonomy corrections and (2) phase space dependent holonomy 
corrections. The effect of inclusion of holonomy correction turns out to be non-trivial and it 
seems that a static solution or a solution with unmodified constraint algebra is inconsistent 
with the inclusion of holonomy corrections in either form. Furthermore, unlike the case of 
inverse triad corrections, here it turns out to be much more difficult to regain the usual 
notion of covariance for the solution of constraints and equations of motion. We point out 
a possible reason for the difficulty and finally conclude in section 5. We will work with 
metric signature ( — h ++) and unless otherwise specified Greek indices {fi, v..) will denote 
spacetime indices taking values (0, 1, 2, 3). 
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2 Classical theory 



The system we are interested in is the minimally coupled Einstein-Maxwell theory with 
spherical symmetry and is described by the action 

S = J d^xiCc + Cem) = J d^x^ (l^^ - J^M^^'^') ' (1) 

where R is the Ricci scalar, F^^, is the electromagnetic field strength, g is the determinant 
of the metric tensor g^^^, and G is Newton's constant. To investigate the effects of loop 
quantum gravity we need the canonical formulation of the above theory which for the 
minimally coupled case can be done independently for the gravitational and the electro- 
magnetic sector. We start with the gravitational sector where, since we are interested in 
loop quantization, we use densitized triads instead of the spatial metric components. For 
spherical symmetry, using the 5u(2) basis Xj this is (see [271 ES, EH] for details): 

E = E^{x)Tssme— + {E\x)n + E\x)t^) smO— + {E\x)t^ - E\x)n) — . 

Momenta conjugate to the triad variables are given in terms of the Ashtekar-Barbero 
connection A^^ = FJ^ + 'jK^, where FJ^ and Kl^ are the components of spin connection 
and extrinsic curvature respectively and 7 is the Barbero-Immirzi parameter [SOj [31] . For 
spherical symmetry 

A = A^{x)T3dx + {Ai{x)Ti + A2{x)T2)de + {Ai{x)t2 - A2{x)Ti) slu + Tg COS M^. (2) 

If one introduces the U(l)-gauge invariant quantities {E'^)^ = (E^)^ + {E'^Y and = 
Al + Al, then the symplectic structure is 

{A,ix),E%y)} = {jK^{x),2E^{y)} 

= Mx),P^{y)} = 2GjSix,y). 

Here the field //(x) is a U(l)-gauge angle and has the conjugate momentum 

P\x) = 2 A^E"^ sin a 

= 4tr ((^Vi + E^T2)iA2Ti - A1T2)) 

(with a defined as the angle between the internal directions of A- and i?-components). 
With this the x-component of the spin connection turns out to be F^, = —t]' implying 
that the corresponding Ashtekar connection is Ax = —f]' + "jK^ (where K^. is the extrinsic 
curvature component). 

In terms of these variables we have the Gauss constraint 

Ggrav[A] = ^y" dxA((E7 + n, (3) 
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the vector constraint 



2G 

and the gravitational part of the Hamiltonian constraint 



(4) 



/^grav[A^] = j dxN\EThK^^E'^ + 2K^K,E- + {l-rl)E'^ + 2r'^E-) (5) 

with = — (£'^)'/2i?'^ the gauge invariant angular component of the spin connection. 
The spherically symmmetric metric in terms of these variables is given by 

ds^ = -N^de + (da; + N'^dtf + E"dl]^ (6) 

Here N and are the lapse function and the shift vector respectively and turn out to 

be the Lagranagc multipliers of the theory while and E^ are the dynamical variables. 
The radial direction is coordinatized by x and dVt^ = d^^ + sin^ 6d(f)'^ is the angular part 
of the metric. In terms of these variables the determinant of the four metric is given by 
= NE^^/E^ sin e. 

We now move on to the Maxwell field for which the Lagrangian density is 

Cem = -\v^gF,,F^\ (7) 

where the electromagnetic field tensor F^y = d^Ay — d^A^ in terms of the electromagnetic 
four potential A^. Momentum conjugate to A^ is 

-'^(-) - 1^ = V^gF^'- (8) 

Since F^" is antisymmetric, we have the well known result 7r°(x) = or in the smeared 
form 

Ci[v] = / dxv{x)T:^{x) w (9) 



which is the primary constraint of the theory. 

Before working out the Hamiltonian we note that in spherically symmetric situation 
there will not be any 9 and components of the electric and magnetic fields. We will 
further assume that there are no magnetic fields. This can be achieved by having the 
Lagrangian depend only on A^ and A^ components of the four potential and with both 
these fields depending only on {t,x). In such a situation only the radial component of tt^ 
will be non-zero 

7r^(i,x) = NE'^Ve^F''^ sine = p^sin^ (10) 
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where, in the above equation, we have introduced the notation = NE'^y/E^F^^. The 
fields Ax and obey the Poisson bracket relation {A^^y) , (z)} = 6^{y,z) or after inte- 
grating over the angular directions we have the relation 



{AM,P%z)} = -S{y,z) 



where y denotes all three spatial coordinates and y denotes only the radial coordinate. 
With TT^(x) = 0, the Hamiltonian 



H 



EM 



d'y 



7r%y)A,{y) - C 



EM 



after integrating over the angular coordinates is given by 

'NE^ip'^f 



EM 



An / dx 



+ P^d,Ao 



_ 2(E^)3/2 

We can do an integration by parts on the last term and write the Hamiltonian as 

'NE^ip^^f 



(12) 



n 



EM 



An / dx 



2{E^f/^ 



Aod^p'' 



(13) 



We already have a primary constraint (|9j) in the theory and following Dirac [52] we 
need to work out the consistency condition(s) that this implies. Since this constraint does 
not involve any other canonical variable apart from = p^sinO, we can forget about the 
gravitational part of the Hamiltonian and calculate the Poisson bracket {p^jHem} which, 
using {Ao{y),p^{z)} = 5{y,z)/ATT, is 



{p'{x),nEM} = {7r°(x),47r j dx 



NE^ip^^f 

_ 2(E^)3/2 



Aod.,p' 



} =p^'(x). 



(14) 



where prime (') denotes derivative with respect to x. For consistency this should equal 
zero, implying that we have a secondary constraint in the theory p^ ^ 0, which in the 
smeared form we write as 

C2[/3] = / dxP{x)p'''{x) ^ 0. (15) 



This would be recognized as the Gauss constraint of the electromagnetic theory. We also 
need to evaluate the Poisson bracket of this secondary constraint with the Hamiltonian to 
see if it gives any further conditions. It is easily checked that this Poisson bracket does not 
give any new condition and turns out to be identically zero (since Hem does not depend on 
A^). We also note that the matter Hamiltonian Hem itself involves the constraint C2[Ao] 
which we separate out and write the electromagnetic Hamiltonian as 



nEM=HEM[N]-C2[A7cAo] 



6 



where 



HEM[N]=A7r j dx 



We now have five constraints in the theory (|3]), Q, ([o]), ( [l5) ) and the Hamiltonian constraint 
H[N] = -^jdxN\E^\-HKlE^ + 2K^K,E'-' + {l-rl)E^ + 2r^E^) 

Note that due to the antisymmetric nature of the Maxwell field and due to the imposition 
of spherical symmetry with Ag = = Afp, there is no matter contribution to the vector con- 
straint. For later purpose (dealing with inverse triad corrections) we also note that in the 
matter part of the Hamiltonian E^ comes with a power of —3/2 while in the gravitational 
part it comes with a power of —1/2. Next we verify that the set of five constraints forms a 
first class system. Before doing that it is useful to note that Ci [u] and C2 [/?] depend only on 
p° and respectively, and in particular are independent of gravitational fields. Similarly 
D[N^] and G[X\ are independent of Maxwell fields. This then immediately implies that 
Ci[i/] commutes with all the constraints. Similarly we find that C2[P] Poisson commutes 
with all the constraints. The gravitational Gauss constraint G[X\ Poisson commutes with 
itself (G[A]) and with the Hamiltonian constraint. Its Poisson bracket with the vector 
constraint is 

{G[X], DIN""]} = -G[A'Ar"]. 

For the vector constraint we have {/^[A^^^], D[M^]} = D[N^M'-^' — N^' M^] and to evaluate 
{D[N''],H[N]} we write H[N] = Hg,^^[N] + Hem[N] and find that 

{D[N%Hem[N]} = Hem[N'N^] + C2[ 



which implies 

{D[N''],H[N]} = H[N''N'] + C2 



(E^)3/2 

AnNN^'E^p''^ 



(E^)3/2 J- 

Finally, we have to check the Poisson bracket of the Hamiltonian constraint with itself - 
{H[N],H[M]}. To evaluate this, as before, we write H[N] = Hg,^^[N] + Hem[N] and 
note that {Hem[N], Hem[M]} = 0, {Hg,,,[N], Hem[M]} = -{Hem[N], H^r.^M]} and 
therefore the only non-trivial contribution comes from the {-ffgrav[-^]) -f^grav[^]} bracket 
and we find 

{H[N],H[M]} = DllE^liE'Py^NM' - N'M)] - G[\E'^\{E'P)~\NM' - N'M)r]']. 

Thus the constraints form a first class system. At this stage we can simplify the analysis 
somewhat by solving the gravitational Gauss constraint G[X\ and the primary constraint 
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Ci[z/]. This removes the pairs (77, P"^) and {Aq, tt") and we are left with the canonical pairs 
{K,{x),E%y)} = {K^{x),2E^{y)} = 2GS{x,y) and {A,{x),p%y)} = ^6{x,y). 

Among the remaining three constraints, while C2[P] and H[N] retain their form (remem- 
bering that now is one of the canonical variables) the vector constraint reduces to the 
diffeomorphism constraint (which we continue to denote by letter D) 



Equations of motion 

The total Hamiltonian Ht = H[N] + D[N^] + C2[P] and using it we can work out the 
equations of motion for the dynamical variables of the theory using Hamilton's equations 
a = {a, Ht}- 



A, 



NK E'f , 

[f^-E^)' + + NK,y/E^ 



/Ex 



, NKl 
N'^K' 



NKlE^ 



N ^ N{E' 



x'\2 



NE^ 



N'E^ 



2(^x)3/2 



+ 



2(E^)2 

N{E'^'y 



NE""" 



+ 



N'WE^ 
+ — ^ + 



E^ ■ 2(E^)3/2 8E^{E^fl^ 2E^^E^ 
N'E""' N'E^'^fE^ 



2{EfY\/E^ Ef 2E'f\fE^ 
f3' NE^p"" 
■4^ + (E^)3/2 



{E^fl^ 







(19) 
(20) 

(21) 



Q^G Z (22) 



(23) 
(24) 



We have six equations of motion and three constraints to solve for six dynamical vari- 
ables and three Lagrange multipliers. Now C2[(3] = implies that is independent of x 
and the last equation above implies that it is also independent of time t thus implying that 
p^ is a constant. We now choose the gauge = and E^ = x"^ and look for static solution 
for the above set of equations. Note that the choice for E^ in conjunction with (6|) implies 
that X refers to the radius of spherical sections. With these choices, (19) and (20) imply 
that for non-zero A^, = = K^^. This then implies that the diffeomorphism constraint 
is already satisfied. We are thus left with the Hamiltonian constraint and three equations 
of motion to solve for A^, E"^ and /3. We will thus have a non-trivial consistency condition 
telling us whether our gauge choice or the assumption of static solution is consistent or 
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not. We find that a consistent solution is given by 



. ,|l_?^ + ^y"'^ (25) 

This determines = Q/\/Au: which we identify as the electric charge and /3 = — \/AixQ / x 
which is the electrostatic potential. Finally substituting the solution (25) in ([6]) we have 

2GM , GQ^\ ^ . , / 2GM , GQ^\-\, , 



ds^ = - 1 + ^ dt^ + 1 + dx' + x'dn\ (27) 



which we recognize as the classical Reissner-Nordstrom black hole solution with mass M 
and charge Q. 



3 Inverse triad corrections 

In this section we work out the consequences of incorporating the inverse triad corrections 
in the Hamiltonian constraint. These corrections arise when loop quantizing the inverse 
powers of triad variables. In the quantum theory it is the flux variables that appear in the 
constraints. These are obtained by integrating the triad variables over two dimensional 
surfaces and their spectrum is discrete containing zero and thus they do not have direct 
inverses. However, by using techniques from the full theory an operator equation for the 
inverse can be defined indirectly which reproduces the classical inverse when quantum 
effects can be neglected [331 El]- We can incorporate this correction at the effective level 
by replacing 1/ ^ a{E^) / E^. The correction function a{E^) can be obtained by making 
use of the techniques referred to above. 

In principle the correction should depend not only on the macroscopic scale as provided 
by E^ but should take into account the underlying discreteness of the theory. This is 
similar to the method of lattice refinement as used in cosmological models [35l [36] . For 
homogeneous or symmetry reduced models the underlying discreteness of the full theory 
is not explicitly visible along the symmetry directions. To remain as close to the full 
theory as possible one then has to include such effects by hand. From the full theory one 
expects that with increasing area size (as determined by E^) the underlying discreteness 
should get refined (since on macroscopic scales the spacetime is smooth). If the underlying 
discreteness is imagined as some kind of a lattice (realized via the spin network states) 
then these considerations would imply that with increasing area the number of discrete 
plaquettes forming a given spherical section should increase. Thus in effect, the number of 
plaquettes is a phase space dependent function which we write as Af{E^). 

In such a situation the correct behavior of the correction function a is determined not 
by the scale E^ but by the scale A = E^ /N'{E^) determining the plaquette size. For 
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A 

Figure 1: The correction function a{A) with A taken relative to A^, := a/7/2£p. 



spherical symmetry we find that the correction function is given by 



VIA + 74/2^V|A- 74/21 
74 



«(A) = 2VA^ ' ' / ' '-^ {2i 



(see [37j for a derivation of the above formula). Fig. [T] shows the behavior of the function 
and we can see that as a; — )■ 0, the classically singular behavior of 1/{E^Y^'^ is regularized 
and in the quantum domain instead of diverging it goes to zero. Thus it is the plaquette 
size which controls the effectiveness of these corrections and for small plaquette size the 
corrections are large. From the figure we also see that the function a(A) — )■ 1 quite fast 
once we are above the scale A^, = (7/2) ^/^£p and the corresponding quantum corrections 
are small. 



We thus correct the classical Hamiltonian (17) to obtain 
H^[N] = J dxA^ 



+ 2ar' - 47rG5E^(p^')^(E^')"^^^ 



(29) 



In the above expression, for generality, different powers of E^ are corrected by different 
a's. We would like to check whether demanding that the constraint algebra be first class 
even after incorporating quantum geometry corrections imposes any restrictions on the 
exact form of the correction functions a and a (relative to a). Note that formally 
can be written as and whereas the first factor will not pick up any quantum 

corrections, the second factor will have quantum corrections due to the inverse power of E^ 
and thus to that extent a is a reflection of quantization ambiguity [3H] . As already alluded 
to while discussing the classical constraints, in the matter part of the Hamiltonian we have 
|^x|-3/2 ^]-^jg Qg^j^not be directly related to lii^^^l"^/^ and in principle the correction a 
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to this factor can be determined independently of a. Thus it would be interesting to see 
whether demanding a first class algebra restricts the function a. 

Since the diffeomorphism constraint does not contain any inverse triad factors and 
because its action is directly represented on quantum states through group averaging, it is 
left unmodified. The Gauss constraint C2\ji\ is independent of triads and therefore it also 
retains its classical form. For this reason the Poisson bracket between any two constraints 
of the pair D\N^\ and C2\P\ has its classical form. Similarly Poisson commutes with 

6*2 since the former is independent of whereas the latter depends only on ]f . We 
therefore only need to evaluate Poisson brackets {D\N^\,}i^\NW and {if*^[A^], iJ'3[M]}. 
After some straight forward but tedious algebra one finds 



Earlier in describing the form of the correction function a(A) we had assumed that 
it depends only on ii^^. A simplistic justification for this assumption can be provided by 
saying that only the inverse triad corrections due to are present and therefore correction 
function should also depend only on this variable. However, in the next section on holonomy 
corrections we will see that such arguments are not always valid. In fact to start with, one 
can allow a to depend on E'p as well. However, as discussed in [13], E"^ dependence of a 
would lead to anomalous terms in the first Poisson bracket above. Demand of anomaly-free 
algebra is thus seen to restrict the form of the correction function. 

Thus we see that for arbitrary a and a we have a first class algebra, the only condition 
being that these be functions of E'^ only. However, even though the algebra is first class, it 
does not have the classical form with the bracket between two Hamiltonians depending on 
a while that between the Hamiltonian and the diffeomorphism constraint depending on a 
both of which are equal to one only in the classical theory. As has already been discussed 
in a previous paper [T^], this implies that gauge transformations are no longer equivalent 
to coordinate transformations. The conventional spacetime picture is thus lost. 

On the other hand we also notice that since the condition of first class algebra does 
not restrict the form of a (and because it is possibly related to quantization ambiguity), 
we can have a = 1, in which case the Poisson bracket between two Hamiltonians has the 
classical form and the usual spacetime concepts hold. Thus if we put the stronger condition 
that not only should the algebra be first class but should retain the classical form then we 
find that a (as also a) is constrained to be unity. Unlike a, however, a is not related to 
quantization ambiguity and thus in general the full constraint algebra would not have the 
classical form. 



With the total Hamiltonian given by = H'^lN] + DIN""] + C2[P], the equations of 



{D[N^],H^[N]} = H^IN'M""] + C2[4'KNN''aE'^p''\E''\-^/% 
{H^ [N] ,H^[M]} = D [a^ {E^KE^^y^ {NM' - N' M)]. 



(30) 



(31) 
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motion are 



NaK E'^ 



/E^ 
aNKl Na 



(32) 
(33) 



2VE^ 2VE^ 8(E^)2V^ 



2{E^y 



2{E^y 



4:{EvyVE^ 



(34) 



{N^'K^y + ,„ - "'Z- ^ - 2NK.„KWE^ 



6a 



2(E--)3/2 

NaE'^ NE^ 6a 



+ 



NaiE'^'f 



iV(E^')2 6a 



^ 6E- 
N'aE""' NaE""" 



2(£'^)3/2 y/E^ 6E'= 8E'/'(E^)3/2 /lE'fiy/E^ SE"" 2Efy/E^ 2E'fiy/E^ 



NaE^'E'^' N"aVE^ N'aE^'VE^ 2iV'E^'v^ 6a 2N E^' E""' ^/E^ 6a 
+ — . + 7T^. + - 



2{E^yVE^ 



Ef 



Ef 



6E^ 



6E^ 



NaiE'-^'f ^ N{E^'y 6a ^ N'aE'^' ^ NaE""" NaE^'E'"' ^ NE"""^ 6a 



+ 



4i?'^(E^)3/2 Efy/E^SE'^ E'p\fE^ Ef\^ {E'PY\f& ' E^ 
Na"y/E^ 



^ ^NaE^ip'^f ^ ^NE^ip^f 6a 
dttG — , „ , ^ h AttG- 



Ev> 



(E^)5/2 



(E^)3/2 6E^ 



p 



13' NaE^p"" 
■4^ + (E^)3/2 







(35) 

(36) 
(37) 



For this set of equations of motion and constraints we now verify that a static solution 
is possible for arbitrary a and a. As before we have six equations of motion and three 
constraints and nine variables to solve for. As for the classical case, the constraint C2W] = 



implies that p^ is independent of x and then (37) implies that p^ is a constant (to be 
identified with the electric charge). Also note that (36), which determines /3 in terms 



of A^, i?*^ and E^ is decoupled from the rest of the equations. For static solution 
is independent of time and its dependence on x is arbitrary which is fine since it just 
determines the electromagnetic gauge. For static gauge we have A^^ = and we choose 
E^ = x^ . Then (32) and (33) imply = = K^- This implies that the diffeomorphism 



constraint is identically satisfied. We are thus left with the Hamiltonian constraint (29) and 



two equations of motion (34) and (35) to solve for two variables A^ and E''^ thus implying 



a non-trivial consistency condition. 

With all this the Hamiltonian constraint (29), after substituting for and its deriva- 
tive, simplifies to 



aE"^ 



X 



ax 2ax 2ax^E'^' 
'E^ ~ ~W ^ ~[E' 



47rG(p^)^aE^ 



0, 



(3^ 
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implying the following differential equation for 

2ax^E'P' + ax^i.E'Pf - ax'^E^ - 2ax^E^ - 4nG{p''fa{E'^f = 0. (39) 



Similarly (34) can be simplified to give a differential equation for 

2ax^N' - ax'^{E'PfN - ax^N + 2ax^N + 2ax^N - A-nG{p''fBi{E'PfN = 0. (40) 
Use of these two equations along with the gauge choice and the solution K^, = = 



shows that the r.h.s. of (35) is identically equal to zero, that is = 0, thus proving that 
a static solution is possible for arbitrary a and a. 

In the next two subsections we will discuss two cases: (I) a = 1 where corrected 
Hamiltonian results in modified dynamics but, with the (spacetime part of) constraint 
algebra being classical, the spacetime properties are unmodified and (II) a = a with a 
correspondingly modified constraint algebra implying that not only the dynamics but also 
the spacetime properties are modified. In this case we will also try and see whether, 
through some suitable modifications, one can salvage conventional spacetime notions. 

At this stage we also have to decide on what form to choose for a. As explained earlier, 
the form of a is determined by using techniques of |33l |3l]. The method works to give 
quantum inverse for (i?^)~^, < p < 1 whereas in the present case we require the inverse 
operator for (£'^)~'^/^ where this method does not work. Since (i?^)~^/^ = ((i?^)~^/^)^, we 
make the obvious choice a = a^. Again there exist different possibilities for the form of a 
depending on how we view the corresponding classical expression, the above choice being 
just one example. There are thus quantum ambiguities present here which will modify the 
exact quantitative results. However, we expect that the qualitative behavior of the solution 
would not depend on such details. 

3.1 Case I: a = 1 

We are interested in static solutions which, with a modified Hamiltonian, will imply quan- 
tum corrections to the Reissner-Nordstrom black holes. At the outset we mention that 
with only one type of quantum gravity correction included, whatever result we obtain 
cannot be considered as giving a complete picture, specially in the deep quantum regime 
where other effects like the holonomy corrections (to be dealt with in the next section) 
or quantum backreaction effects (not studied here) might become dominant. The present 
analysis should be seen as an exploration of the kind of quantum effects that result from 
different corrections and which correction(s) dominate in different regimes. At the present 
level of development, where not much is known about the quantum behavior in the full 
theory such models can provide intuition for further developments. 

With the gauge choice = and E^ = x^ already made we need to solve equations 



(39), (40) and (36) to determine ii^*^, N and /3. Furthermore, as mentioned in the beginning 
of the section, we need to include the effects of lattice refinement in our solutions such that 
the correction function a = a{E^,J\f). In this sub-section we will consider two cases - 
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(i) M = const corresponding to constant patch number and (ii) M (x x"^ corresponding to 
non-constant patch number. 



With a = 1 and a = a^, equations (39) and ( [40| ) become 

2x^E^' -2x^E^ -ax^E^ + ax\E^f -47iGa^{p'')\E^f = 0, (41) 
2x^N' - ax^iE'^fN - ax^N + 2x^N + AnGa%p''y{E'^yN = 0. (42) 

To solve these equations we first note that classically E"^ = xj {l — 2GM/x + GQ'^/x'^y/'^ 
and = Q/\/4Tr. Also for this class of correction, the modified Schwarzschild solution as 
found in |T3] gives Ef = x/{l- 2GMf^{x)/xY'^ where the function /„ (x) is a solution 
of the equation xf'^ = (1 — a)/^ with the subscript in signifying the dependence of 
the functional form of / on the refinement scheme. Since we know that for Q = our 
solution for E'^ should go over to that for modified Schwarzschild and that in the absence 
of quantum corrections it should go over to the classical solution, we make the ansatz 
E'f = X / {1 — 2GM fa{x) / x+GQ%a{x) / x'^Y^'^ and substitute it in (41) to obtain an equation 
for ha{x): 

xb'^ - 2ba + aba + = 0. (43) 

Before solving this equation we find the form of the lapse function N as well. We again 
make an ansatz which is motivated classically and also from the inverse triad corrected 
Schwarzschild solution and write N = d„(x)(l - 2GMfa{x)/x + GQ%a{x) / x^f/'^ / fa{x) 



with da{x) to be solved for. When used in (42) one finds da{x) to be a constant which we 



choose to be unity to get the correct classical limit. The metric for this guage is then given 
by 

\ X X^ J \ X x^ J 

In the above solution we note that function fa multiplies the mass M and function ba 
multiplies the charge Q. Using the language of perturbative quantum field theory, these 
can be thought of as the mass renormalization and charge renormalization, respectively, 
due to quantum gravitational effects |^ The metric coefficient gu is also corrected by 
and is akin to wave function renormalization. We further note that at the present effective 
level, Newton's constant G does not seem to be renormalized due to quantum gravity 



effects. Using the expressions for E''^ and in (36) we get the electrostatic potential 

15 = v^Q / dx-^ , (45) 

J JaX 

and the factor a^/fa contributes to the renormalization of the electrostatic potential. It 
is interesting to note that the renormalization of the electrostatic potential has different 
factor compared to that for electric charge. 

With only the Hamiltonian constraint modified and the constraint algebra retaining 
its classical form, we expect that solutions of equations of motion and constraints will be 



^Terminology suggested by Roniesh Kaul when shown some of the solutions. 
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mapped to other solutions under general coordinate transformation. As in [15], we verify 
this by constructing the Painleve-Gullstrand analog of the above metric. Since the metric 



in (44) is time independent, = (1, 0, 0, 0) is a Killing vector and therefore we can choose 



QapW^it — ~1) where = da;"/dT is the four- velocity of the observer falling freely from 



infinity (starting at rest) in the Reissner- Nordstrom metric (44). T is the proper time 
along this trajectory and serves as the time coordinate for the Painleve-Gullstrand metric. 
We also have Qapu'^u^ = — 1- Using these two equations we can solve for which, with 
dT = — Uadx", leads to 

iT^atJl-^-^^^Y (fl-l^^^-^)"\.. (46) 



This when solved for dt and used in (44) gives the Painleve-Gullstrand metric 

= -dT^ + I dx + V - 1 + - ^5^dr^ + ^'d^]l (47) 

Comparing this with (g we have = 1, A^^ = (/^ - 1 + IGMfJx - GQ^^jx^fl^, 
E"^ = x/ fa and = in the Painleve-Gullstrand gauge. It is straight forward to check 
that this solution satisfies all the constraints and the equations of motion. This once again 
confirms that for unmodified constraint algebra even if the dynamics is modified due to 
the modified Hamiltonian, solutions of constraints are mapped to other solutions under 
general coordinate transformations. 

Having verfied that the spacetime picture is valid in the present case we can now analyze 



the metric in (44) for horizon properties. We find that as opposed to the classical case, the 
metric coefficients gu 7^ g'^l- However since in the factor > (it being zero only at 
X = 0), the location of the horizon is given by the solution of 

1-^ + ^=0. (48) 

X x^ 

To explore the horizon properties we will need to consider the explicit form of functions 
foL and ha for different refinement schemes. 

Constant patch number: For M = constant, the solution for fa as found in [15] is 



Ux) = ^jj^ , x'>Afjil/2m 



1 /2 

X + v/^^^Ar^472) [x + y/x^+M^£l/2^ " 

9p-7r/4™„{l-Q(z))/2^| arctan(^ x2/(A^7<??./2-x2)) 

fa{x) = , , x^<Ar^il/2 (50) 

(Ar74/2)V4(a; + + AA74/2)i/2 ' 



which is plotted in Fig. [2j Equation (|43j) is difficult to solve analytically and we will 
therefore work with its numerical solution. For the purpose of numerics we choose Af = 1 
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Figure 2: Function fa{x] 
tive to X* := ^/'yAf/2ip. 



with X taken rela- 



Figure 3: 
tive to X* 



Function ba 



X 



with X taken rela- 



setting G = 1 and measure x in the units of (7£p/2)^/^. Since in the classical limit b{x) — )■ 1 
and because the function a (A) — )■ 1 very fast for A > A^,, we choose the boundary condition 
6(100) = 1. The behavior of b{x) is shown in Fig. |3} To analyze the horizon condition, we 
rewrite (48) as (with G = 1) 



M 



X' 



2xfa 



(51) 



For a given value of M and Q the above equation will determine the location of the horizon. 
Since we do not have analytic solution for b^, we evaluate the horizon condition graphically. 
If we plot the r.h.s. of ( [sT] ) (which we call the horizon curve) as a function of x, then the 
points where the line M = constant intersects this curve gives the location of the horizon. 
This is shown in Fig. |4]for Q = 2, where for comparison we have also plotted the classical 
horizon condition M = (x^ + Q'^)/2x (shown dashed in the figure). 

We note that though the classical curve has only one extremum corresponding to ex- 
tremal black holes, for the quantum corrected case we have two extremas. To reduce the 



dependence on the results of numerics, we extract more information directly from (51) 
by taking the x-derivative of the r.h.s. which works out to a(x^ — a'^Q'^)/2fax'^. Since 
a/2fax'^ > 0, we can extract the desired information from the factor (x^ — a^Q^), which 
is plotted in Fig. |5] and this again shows that there are two extremums. One of these 
extremums corresponds to the classical value while the other is in deep quantum regime. 

Interestingly the presence of non-classical extrema leads to the possibility of the black 
hole having three horizons as can be seen from Fig. |5] where we find that the slope of the 
horizon curve on the left of the inner (non-classical) extremum is positive implying that 
for the choice Q = 2, there are values of the mass of the black hole which can lead to three 
horizons (this feature is difficult to make out from Fig. |4] since the slope of the curve in 
this region is very small). 

Another very interesting feature is that there now seems to be only a finite range of 
mass, approximately 2 < M < 3, for which there exist two (or possibly three) horizons. 
For M > 3 there exists only one horizon (corresponding to the outer horizon of classical 
black hole) and thus one of the distinguishing features of the Reissner-Nordstrom black 
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Figure 4: The r.h.s. of (51) (called the hori- 
zon curve in the text) for Q = 2 with x taken 
relative to := a/7A/'/2£p. 



Figure 5: The factor (x^ — a^Q^) (in the 
derivative of r.h.s. of ([sI])) vs x for Q = 2 
with X taken relative to x* := \/'^M jU^. 



hole is lost. Since in deep quantum regime other effects, like holonomy corrections and 
quantum backreaction effects which are not considered here, can play a significant role we 
do not analyze the implications of the presence of the extra horizon (and for the same 
reason do not attempt to draw a Penrose diagram corresponding to this solution). 

To explore the dependence of the horizon properties on the value of the charge of the 
black hole we now consider the case where the charge of the black hole is small, say Q = 0.1. 
The relevant plots in Figs. [6] and [7] now show different features and we find that for inverse 
triad corrected solution there is now only one horizon for any value of M > 0.34 (in Planck 
units with 7 absorbed) and no horizon forms below this value. This is in agreement with 
the results in [39l [15] where there is a mass threshold for the formation of horizon. 

Furthermore there is no extremal solution possible in this case. For the classical black 
hole on the other hand, the extremal solution corresponds to M = 0.1 and for mass greater 
than this value there are two horizons. For smaller values of charge Q, one is tending to 
the Schwarzschild limit of the Reissner-Nordstrom black hole and the mass threshold found 
above is comparable to the mass threshold for the formation of Schwarzschild black hole 
with inverse triad corrections |15j . 

To find the largest value of Q for which these features persist we note that for non- 
existence of extremal solution, the condition (x^ — a^Q^) > should hold for all x. For 
given form of a this implies that for Q < l/\/2 we have features similar to Q = 0.1 while 
for Q > I/V2 one gets the features discussed above for the case Q = 2 except that for 
large enough values of Q such that (x^ — a^Q^) < the inner extremum disappears and 
the slope to the left of the (outer) minimum in Fig. |4] remains negative and therefore we 
do not get three horizons. However, there is still a mass threshold beyond which the inner 
horizon disappears. 

In all the cases discussed above, one notes that the geometry for the quantum corrected 
case quickly merges with the classical geometry as one moves beyond the Planck scale and 
that all the interesting features correspond to deep quantum regime where, as already 
mentioned, role of other quantum effects needs to be taken into account to draw reliable 



17 



Q=0.1 









Q=0.1 /' 























Figure 6: The r.h.s. of (51) (called the hori- 
zon curve in the text) for Q = 0.1 with x 
taken relative to := a/7A/'/2£p. 



Figure 7: The factor [x^ — a^Q"^) (in the 
derivative of r.h.s. of ([5l|) vs x for Q = 0.1 
with X taken relative to := a/7A/'/2£p. 



conclusions. We also find that the singularity is not resolved with Ricci scalar diverging at 
X = (classically, Ricci scalar is zero and it is the Kretschmann scalar which diverges). 

Non-constant patch number: For non-constant patch number we choose a power-law 
ansatz such that Af oc x^. This includes the special case where the patch size A = /M 
is a constant corresponding to p = 2. In general we can write p = 2 + e with e = 
corresponding to the case with constant patch size. The correction function a in (28) takes 
the form 

J\l + aWx'\- J l-a?h'^x'\ , , 

" = ' ('2) 

where a} = 7^p/2 and 6 is a dimensionful constant needed to make A/" = iP'x'^^'^ dimen- 
sionless. As has been observed in p!S], only for e < do we have physically reasonable 
solutions. For e > 0, the patch size becomes smaller as we move to orbits with larger 
radius (and the correction function a tends to zero for large x instead of tending to one), 
implying that quantum effects due to inverse triad corrections become dominant which is 
unphysical. We now have three parameters in the problem 6, e and Q and numerically it 
is difficult to explore the full parameter space. We looked at a few representative cases on 
which we now briefly comment. 

For h = 0.2 (a = 1 in the chosen units), e = — 1 and Q = 2 the behavior of the 
solution is very similar to the quantum corrected case with constant patch number with 
Q = 2. Similar behavior is observed when 6=1 other parameters remaining the same. 
Also for Q = 0.1 the inner horizon seems to disappear which was also the case for constant 
patch number. This is consistent behavior since in this case the patch size is growing with 
increasing x as is also the case when the patch number A^ is constant. On the other hand 
when e is close to zero, say e = —0.1, the behavior is very similar to classical solution for 
different values of b and Q. 

When e = 0, then the correction function a is a constant and the requirement that in 
the classical regime it should be close to one implies that ab <^ 1 and it turns out that a 
is close to but greater than one. In this case we can analytically solve (43) for ba (where 
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we do not substitute the form (52) but use a itself remembering that it is a constant): 

= cix^-" -, (53) 

a — I 

where Ci is a constant. Since a is close to (and greater than) one, we find that imposing 
the desired hmit 6^ — )■ 1 for x — )■ oo forces one to choose Ci = = 6. For this case the 
function /„ is [I5] 

/„(x) = C3a;^-^ (54) 

and we see that is very slowly decaying function of x remaining almost constant for 
a large range of radial coordinate x. Thus, in general, both fa and ha do not have the 
expected large x behavior of tending to unity. Both the functions however are very slowly 
varying functions of x (for suitably chosen ci and C3) and as for the quantum corrected 
Schwarzschild in [15] lead to a metric which in the asymptotic limit is conformally fiat. 



3.2 Case II: 



a = a 



We now put a = a everywhere in the set (39) and (40) (remember that already the gauge 
choice = and = x^ has been made in writing these equations). Proceeding as 
in the previous case we make the classically motivated ansatz E"^ = x/{l — 2GM/x + 
GQ"^ ga{x) / x'^Y^'^ (for this version of inverse triad corrections it turns out that there in no 
correction function in the mass term, see [15]). Substituting this in (39) we get 



Xga = 9a 



a^. (55) 
To solve for we make the ansatz N = ha{x){l — 2GM/x + GQ'^ga{x) /x'^Y^'^ which on 



substituting in (40) and making use of l\55u gives 



ha 



a 
a 



(56) 



with the solution ha{x) = 1/a (where the constant of integration has been chosen to be 
1). With this it turns out that 



13 = V^ttQ / dx 



a 



X 



2 ■ 



(57) 



Compared to the classical solution there is an extra factor of in the integrand which can 
be thought of as the renormalization of the electrostatic potential due to quantum gravity 
effects. 

With the phase space solution known, the next task would be to write a spacetime 
metric corresponding to this solution. The obvious choice for the metric based on the 
solution of phase space variables would be 



^s' = - 



1 



a' 



1 - 



2GM GQ^go 



X 



X 



x" 



dt^ + 1 



2GM ^ GQ'gajx] 



X 



dx^ + xMQ\ 
{5i 
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However, as has aheady been noticed in [15] in the context of vacuum Schwarzschild 
solution as well as for perturbative matter contributions, such a solution in not covariant 
(the slash in /Is^ indicates that this is a formal construct and does not correspond to 
a covariant object). By this we mean that proceeding analogously to the previous sub- 
section if we perform a coordinate transformation on this metric to go over to the Painleve- 
Gullstrand like coordinates we find that the new metric does not satisfy the complete set 
of equations of motion and constraints. Specifically, following the standard procedure one 



finds that in terms of the Painleve-Gullstrand like coordinates the metric in (58) becomes 



from which comparing with (g we identify N = I, = x/a, N"" = {a^ - I + 2GM/x 




X 



2 



GQ^g^{x)/x^y/^ and 

Even though this form for the phase space variables (as also the metric in (59)) has 
the correct classical limit, these functions fail to satisfy the Hamiltonian constraint. Thus 
we see that for the case of modified constraint algebra, solution(s) of constraints (and of 
equations of motion) do not map to other solutions of constraints under general coordinate 
transformation. This has to be contrasted with the previous case for a = 1, where the 
classical constraint algebra applied and hence solutions of constraints mapped to other 
solutions under general coordinate transformation. 

With quantum corrections modifying the constraint algebra, we therefore seem to loose 
one of the most crucial aspects of the classical theory where gauge transformations in phase 
space correspond to coordinate transformations in spacetime and solutions of constraints 
are mapped to other solutions. In such a situation where conventional spacetime notions 
loose their meaning one cannot define concepts like metric, black holes and their horizons, 
which are used all the time in the classical theory. 

However, we also note that the constraint algebra is only deformed and has not become 
anomalous. Thus we expect that some notion of covariance should still survive. It is a 
non-trivial task to identify what this alternate notion could be, though whatever it be it 
certainly has to go over to the classical notion in the absence of quantum corrections. Since 
in the quantum theory a priori there is no spacetime but only phase space, it is possible 
that the mapping between the phase space and the spacetime that one is familiar with from 
the classical theory might only be a classical realization of a more general mapping. Below 
we provide one such alternative. However, at the outset we mention that the suggested 
mapping has not been obtained in a systematic manner and it is not clear how general it 
is. On the other hand at the present stage of developments it is also useful to know various 
examples where similar methods work since they can provide hint for future developments. 

For the rest of this section we will need to distinguish between the lapse function 
occuring as part of the phase space A^ps (obtained from solving the constraints and equations 
of motion) from the lapse as it is used to describe the metric A'^met- We will show that 
although classically A^ps = A^mct, for the case with inverse triad corrections where a = a, a 
modified mapping N^^t = «A^ps leads to a covariant metric. We have the following solution 
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for the phase space variables: 



2GM GQ^g^{x)Y^'^ _^ f ^ 2GM GQ^g^{x^^^'^ 



E^ = x[l + " 7^ ' , iVp, = a-Ml + , . 

\ X X'^ J \ X X^ J 

(60) 

This would imply that under new mapping A^met is given by 



X x^ 



^^^^ 2GM^«?^ . ,61) 



With = x^ and A^^ = we then have the 'metric' corresponding to the classical 
Schwarzschild-like coordinates 



d.^ = - ( 1 - — + dt^ + f 1 - + I d." + xMff . (62) 

X x^ I \ X x^ 



Next we construct the Painleve-GuUstrand analog of the above metric. Following the 
procedure of the previous sub-section we find 

dT = d«+fl-H^ + 5«%MV7H^-^«%MV^x. (63) 

\ X x^ J \ X x^ J 



which when solved for dt and used in (62) gives the Painleve-GuUstrand metric: 

2 



/2GM GQ^gaix)^^^^ 
\ X x'^ 



ds^ = -dT^ +\dx+{ - ^ 7^ ^ dT + x^dn\ (64) 



From this we get N^^g = E'^ = x, = x^ and A^^ = {2GM/x - GQ'^ga{x)/x^y/^ . To 
go to the phase space we again need to use A^met = ^A^ps, which implies that A^ps*^ = a^^. 
It is now straight forward to check that this choice for the lapse function along with the 



form for E^, E''^ and A^"^' when used in equations (32) and (33) determines = —N^ 
and = -A^^ . It can now be verified that this solution satisfies the diffeomorphism 
constraint and the Hamiltonian constraint as well as the remaining equations of motion 



(34) and (35). 



It is interesting to note that putting Q = in (62) or (64) we get back the classical 
solution exactly and there is no trace of inverse triad corrections. Thus if this modified 
mapping is used for vacuum solutions one would think that the covariance is recovered 
because one is effectively ignoring the quantum corrections using the new mapping. With 
Q 7^ we clearly see that this is not the case since now there is non-trivial quantum 
correction present in the form of function ga{x)- 

Here we therefore have an example where, even though the constraint algebra is modi- 
fied, one could still find a suitable notion of covariance. As required, the new mapping goes 
over to the classical mapping in the absence of quantum corrections since then a — )■ 1. In 
the next section on holonomy corrections, where again the constraint algebra gets deformed. 
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Figure 8: Fun ction g ^jx) with x taken rela- figure 9: The r.h.s. of 1^ (called the hori- 
tive to X ■= \p^^(Wl2l-p curve in the text) for Q = 2 with x taken 

relative to x* := a/7A/'/2£p. 

we will however see that in general it is much more difficult to identify an alternative notion 
of covariance and that probably we are lucky for the case of inverse triad corrections. 
With a covariant solution in hand, it is now meaningful to look at the properties of 



the solution. From (62) we find that compared to the previous case where a = 1, in the 
present version of inverse triad corrections (and for the new mapping between phase space 
and spacetime) only the electric charge is renormalized and that there is no mass or wave 
function renormalization. However, as in that case, the gravitational constant is still not 
renormalized. We now look at the horizon properties for this version of the inverse triad 



corrections also. From (62) we see that the horizon is given by the solution of 



l-?^ + ^%M.o, (65) 
which, as in the previous case we write as 

M=^^i^. (66) 

2x ^ ' 

The behavior of function is shown in Fig. [8] and we see that it is very similar to and 

To analyze the horizon properties we proceed analogously to the previous case. We 



note that equation (66) is similar to (51) except that there is no /q, in the denominator. 
To explore the features of the horizon, as before we take derivative of the r.h.s. of the 
above equation with respect to x which evaluates to [x?' — a^Q"^) /2x^ . This is exactly the 
expression used to discuss the properties of the horizon in the previous case (and shown 
in Fig. |5| and thus the conclusions drawn for the previous case where a = 1 continue to 



hold in the present case also (though the exact values will differ). For this reason we do 



not repeat the detailed analysis here except for showing a plot for the r.h.s. of (66) for the 
case Q = 2 in Fig. |9| 

In this figure the inner extremum and the corresponding non-classical region to the left 
of it can be seen very clearly. We can also explicitly see the presence of the three horizons 
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for a finite range of mass. As an example we have shown one such case where the horizontal 
line corresponding to M = 2.8 intersects the curve at three different points. Again we note 
that the inner horizon disappears for mass greater than some finite value (M ^ 3.2 in the 
case under consideration). Also for this version of inverse triad correction, the singularity 
is not resolved, though (un)like the (previous) classical case the Ricci scalar continues to 
be zero and it is the Kretschmann scalar which diverges. 

4 Holonomy effects 

We now move on to consider the effect of including holonomy corrections. These arise 
because in the quantum theory there is no representation for connection components as 
operators but only their holonomies have well defined operator analogs. Compared to 
inverse triad corrections these corrections are of a different nature and are similar to (but 
not the same as) the higher derivative (higher curvature) terms in the effective theory. 
These corrections have not been analyzed much for spherically symmetric backgrounds 
though this has partly been done for marginally bound Lemaitre-Tolman-Bondi models in 

In the spherically symmetric background there are two independent curvature compo- 
nents - K^p along the homogeneous angular directions and along the inhomogeneous 
radial direction. At an effective level that is being attempted here, one can incorporate 
holonomy effects by making a replacement — )■ sm{6Ktp)/6 in the Hamiltonian, where 6 
refers to the path length used to calculate the holonomy and is a reflection of the scale of 
discreteness of the underlying quantum geometry. Since corresponds to homogeneous 
directions, the holonomies corresponding to this component are point holonomies familiar 
from loop quantum cosmology. 

This form of the function takes into account the fact that point holonomies correspond 
to replacing connection components by periodic functions in the Hamiltonian (the exact 
form of which will be found when fully working with the quantum theory rather then 
the effective approach being followed here) as also the fact that in the classical limit 
where the discreteness scale 5 — )■ 0, we should recover the correct classical expression 
sm{5K^)/5 K^. 

Incorporating the effect of holonomies is more non-trivial since they refer to the 
inhomogeneous radial direction. At the effective level one can again attempt a replacement 

— )■ sin(/i^'^)// where / is the path length along the radial direction. However because 
the radial direction is inhomogeneous, unlike the scale 5 used for K^p holonomy the scale 
/ = /(x), there being a non-trivial dependence on the radial coordinate which needs to be 
handled carefully to ensure that the diffeomorphism invariance of the theory is not lost. 
Below we will only include the effects of the K^^ holonomy leaving the more complicated 
case of the holonomy for future work. 

As already mentioned, consistent inclusion of K^^ holonomies could be achieved in the 
marginal LTB model. However a more systematic treatment for their inclusion has been 
given in [25] working directly with constraints and obtaining the form of corrections by 
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demanding that the first class nature of constraint algebra is preserved even after including 
these corrections. We will follow the same approach in the following. 



4.1 Phase space independent holonomy corrections 

As before we have three constraints in the theory and we will correct only the Hamiltonian 
constraint (17), repeated here for convenience 

H[N] = -^JdxN\E'=\-^KlE'^ + 2K^K,E^ + {l-Tl)E^ + 2T'^E'^) 

(67) 

We see that occurs at two places in the above expression and that these occurences 
are with different powers in the exponent. A priori there is no reason that the replacement 
K^p — )■ sin(5i^'(p)/5 at both the places will give a consistent picture. Thus to begin with 
there is a two-fold functional degree of arbitrariness. To be general, therefore, we do not 
fix the form of the correction functions from the beginning and replace the quadratic K"^ 
in the first term by ff{Kip) and the linear K^p in the second term by f2{K^p) (for the ease 
of notation we will not show the explicit depedence of /i, /2 on K^p in the following). With 
this the Hamiltonian constraint becomes 

H'^m = j dxN\EX'^{flE'P + 2f2K^E^ + {l-Tl)E'^ + 2T'^E^) 

As for the case of inverse triad corrections, with only the Hamiltonian constraint re- 
ceiving corrections, we only need to evaluate the Poisson brackets {D[N^],H^[NW and 
{H^[N]^ H'^[MW. Again a lengthy exercise in algebra gives 

{D[N% H^[N]} = H^[N'N'^] + C2[47rArArV^^(^^)"^^1, (69) 

which is the same as the classical expression. The other bracket gives a more non-trivial 
result 

{HQ[N],HQ[M]} = D[^E^{En-'{NM'-N'M)] + ^ J dz{NM'-N'M)^ (^f, - f,^ 

(70) 

Clearly for a first class algebra the second term on the r.h.s. should vanish. This happens 
when the integrand in that term vanishes 

This gives a non-trivial condition on the form of the functions thus confirming that one 
cannot make the replacement K^p — )■ sin{6K^) / K^p everywhere. To solve the above equation 
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we can choose a form for either fi or for /2 and determine the form of the other function. 
Thus at the effective level there are several choices. However, it is interesting to note that 
demanding a first class algebra restricts the arbitrariness to just one functional degree 
of freedom. This we can fix by appealing to the full theory where point holonomies will 
appear as periodic functions. Since /2 appears as a replacement for the linear in the 
Hamiltonian, it seems more appropriate to choose f2{K^) = sm{6K^)/6 and determine 
the form of /i rather then the other way round. With this choice for /2 we find that the 
algebra remains first class if fi{K^) = 2sm{6K^/2)/6. In general, /2 = sm{nSK^)/n6 
implies /i = 2sm{n6K^/2)/n6. 

Thus we find that like with inverse triad corrections a first class constraint algebra 
is possible even with holonomy corrections. However, unlike the case of inverse triad 
corrections, we note that in the present case we cannot obtain an unmodified algebra. To 
obtain an unmodified algebra we would need to impose the condition df2/dK^ = 1 (see 
(70)) which imphes /2 = and (71) then leads to /i = K^. Retaining classical algebra 
thus forces one to give up the holonomy corrections. 

In the full theory where all the corrections would be present simultaneously (including 
the holonomy corrections), one would therefore expect that even if one could restrict the 
constraint algebra to be first class, it would be highly unlikely that it will retain its classical 
form. And it has already been observed that with a modified constraint algebra there is 
no standard spacetime interpretation of the solutions to the constraints and the equations 
of motion. 



4.2 Phase space dependent holonomy corrections 

In the previous section we assumed that the holonomy corrections are independent of phase 
space. In principle, just like inverse triad corrections, these corrections should also be phase 
space dependent [35l |36]. We have seen that 6 refers to the scale of discreteness in the 
quantum theory. At the classical level one would expect that this discreteness should tend 
to zero relative to a classical length scale. For spherically symmetric backgrounds under 
consideration, provides one such length scale, referring as it does to the area of spherical 
spatial sections for each choice of x. Thus one expects that the discreteness scale 6, instead 
of being a constant should be a function of so that with increasing E^, 6{E^) tends 
towards zero. 

From the perspective of the full theory of quantum gravity, where the underlying quan- 
tum geometry is represented in terms of spin networks, this corresponds to the fact that 
under a Hamiltonian evolution new edges and vertices are created (and also destroyed) 
thus refining the graph underlying a given state of the theory. In loop quantum cosmology 
the fi scheme |10] provides one realization of this effect. (One can get phase space depen- 
dent holonomy corrections even without lattice refinement if one combines them with the 
inverse triad corrections and demands a first class constraint algebra PS].) 

For this section we therefore have /i = fi{K^,E^) and /2 = f2{K^, E^). We again 
demand that with this choice of correction functions the constraint algebra be first class. 
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Only the {H'^[N], H^[M]} is non-trivial and gives 



{H^[N],H^[M]} 



D[^E''{E^)-'^{NM' - N'M)] 

(JJ\ 



+ 



2G 



dz(NM' - N'M)^ (f2-fi^ + 2E'^ 



.(72) 



For a first class algebra we therefore require 



(73) 



and compared to (71) there is an additional term. To solve this equation we choose 
f2{K^p,E^) = sm{6{E^)Kip)/6{E^). As before, a basic condition on such correction func- 
tions is that in the classical limit /2 — )■ K^^. Following the arguments in the beginning of 
this sub-section and noting that since is dimensionless so that 6{E^) is also dimension- 
less , we choose 6{E^) ~ {•yip/E^y, p > 0. This choice automatically gives us the classical 
behavior for E^ ^ 'yip. With this choice for 6{E^) we have 



E^ 



sm 



E^ 



(74) 



and from the first term of (72) we note that with a non-trivial K^p dependence of /2, 



undeformed constraint algebra is not possible. This form for /2 when used in (73) gives 



two solutions (differing in an overall sign) of which the physically important one is 



(2) 



1/2 



E^ 



+Cp{E^) 



1/2 



1 + 4p) 



cos 



74 
E^ 



+ 2pKp sin 



02 \ V 

'P 



(75) 



where in the above equation Cp{E^) is an arbitrary function (note that (73) involves partial 
derivative of /i only with respect to K^p and therefore does not completely fix the E^ 
dependence). We can fix the form of Cp{x) by demanding that /i — )■ when yip/ E^ — )■ 0. 
Doing a reduction consistently at order i^"^, we find that Cp{E^) = (1 + 4p){E^yp / {yipy^ 
gives the correct classical limit. 

In the following we will mostly work with an arbitrary value of p. However certain 
well motivated considerations lead one to choose specific value (s) for p. In the previous 
section on inverse triad corrections, the effects of lattice refinement were encoded in the 
number of patches J\f{E^) which was seen to be related to the area of spherical sections. 
Similarly in the present case we know that 6{E^) refers to the length of the curve on 
the spherically symmetric orbits along which the holonomy is evaluated. The lattice gets 
refined as one moves to spherical sections of larger radius and one possibility is that the 
number of lattice sites on a given orbit increase in proportion to the radius (corresponding 
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to constant lattice size) and thus 6{E^) itself should vary inversely as the radius implying 
that we choose p = 1/2. For this value oi p (75) gives 



1,1/2 



(2) 



1/2 



7^ 



cos 




f E^ 



sm 




1 1/2 



(76) 

(even though 



At this point it is useful to note that ( 75 ) and ( 76 ) are not periodic in 
f2 is) unless p = 0. It is therefore not clear how such corrections can arise from the 
underlying theory where it is expected that point holonomies should appear as periodic 
functions |^ On the other hand putting p = in the equations would kill the phase space 
dependence of holonomy corrections. Demand of a first class constraint algebra for phase 
space dependent holonomy corrections thus seem to imply non-trivial corrections from the 
perspective of the full theory. 

We now look at the dynamics for this class of corrections. The Hamiltonian constraint 



is given by (68) (but we remember that fi and /2 in that expression now depend on both 



and E^) and the equations of motion are 



A, 



N^E"" +2Nf2VE^ 



{N'^'E^py + 



NhE^ df 



Nfl N 



- + NK^^E^ 



NiE'-'f N'E'^'^/E^ ^Nip^^f 
H : . „ , — H xn 1- 2itG- 



2^f& %{E^Y^fE^ 2{E- 



(^x)3/2 



(77) 
(78) 

(79) 



Nf^E'P 2NfiE^ dfi NfiK^ ^^^^ rF^df2 NE^ 



2(E^)3/2 

iV(E^')2 ^ NE'^ 



NE^'E""' 



'E^ 

N"y/E^ 
+ — ^ + 



dE^ 2(E^)3/2 
N'E""' 



8£'¥'(£^^)3/2 2E'^\/E^ 2{Efy\^ E'p 2E'^^/E^ 



N'E^ JE^ 



V\2 



6nG 



NE^ip'^f 
{E'^f/^ 



(3' NE^p"" 

+ (^^•)3/2 





^0) 
^2) 



As in earlier cases, we try to solve the constraints and the equations of motion looking 



for static solutions with A^^ 



X 



and as mentioned earlier, we do not fix a value for 



the refinement parameter p. However, in the present case we find that for this gauge choice 

= 0. This implies that for non-zero A^, /2 
n G I. In such a situation even /i 



(77) gives 2Ar/2. 



K^p = niTX 



^7(74 



E^ 

2p 



which is possible only if 
and except for the non-classical 



form of Kip for n 7^ 0, nowhere else in the solution do we see the quantum effects. This is 
so because putting /i = = /2 everywhere in the constraints and in equations of motion 



^Author is thankful to Martin Bojowald for pointing this out. 
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is equivalent to going back to the classical solution with K^^ = = K^- This happens also 
for the case where we consider phase space independent holonomy corrections. 

We therefore find that choosing A^^ = and demanding a static solution always gives the 
classical solution unless there is some non-trivial correction involving in the Hamiltonian 
(since (77) involves derivative of the Hamiltonian with respect to Kx)- Thus K^p holonomy 



corrections alone seem to preclude a static solution as well as an unmodified version of 
constraint algebra. However a stationary solution might still be possible, which we look 
for next. 

To obtain time independent solution, we choose = x"^ but keep the functional 
dependence of other phase space variables arbitrary. As before, C2[f3] = and = 



together imply that is a constant which we again choose to be Q/y/Air. Also (81) 
determines /3 in terms of E^, and p^. Equation (77) implies A^^ = —Nf2- Solving the 
diffeomorphism constraint gives = K'^E^'^/x. Using these in (78) then gives the relation 
E''^ = cx/N, c a constant. This is the (classically) expected relation between E''^ and A^ 
(eg. the classical Reissner-Nordstrom solution) and we will choose c = 1. Using these 



results in (79) and in the Hamiltonian constraint (68) shows that these two equations are 



equivalent. Thus we are left with (80), (81) and the Hamiltonian constraint to solve for 
K^p, N and /3. The Hamiltonian constraint gives 



fl 

X 



+ V2K 



1 

X 



?>x 



(E 



ip\3 



^3) 



Replacing E'^ hj x/N this equation can be written as 



X 



1 

2hK'^ = -- + — + 2NN' 



GQ' 



x-^ 



^4) 



We note that the left hand side of this equation depends only on K^, and the right side 



only on A^. Next we write (80) entirely in terms of K^p and A^ to get 



2x2 



-2/1 



dh hK 



dE^ 



-2xK' 



X 



dE' 



+ 



Ar2 NN' 



2x2 2x2 



X 



-NN"-{N'y^ 



3GQ^ 



2x4 
(85) 

As for the Hamiltonian constraint, this equation separates into two parts - left hand side 
depending on and the right side depending on N. It turns out that the right side of 



the above equation is nothing but one half of the x-derivative of the r.h.s. of (84) and can 



therefore be replaced by one half of the x-derivative of the l.h.s. of (84) to get the equation 
entirely in terms of K^^. Finally substituting for fi and /2 from (75) and (74) along with 



the required derivatives one finds that the above equation is identically equal to zero, that 
is, Kx = is identically satisfied. 

Thus with the choice E^ = x^ and the assumption of time independence we have solved 



all the equations except for (84) with A^ and still undetermined (of course (81) is also 
there but it just determines /3). We therefore have the freedom to arbitrarily choose one 
of the functions. We now present the solution to (84) which classically would correspond 
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to the Reissner-Nordstrom solution in (27). To obtain such a solution we make the choice 



2GM GQ'^^^'^ 



N=[l + ^ . (86) 



This using E'^ = x/N implies 



E'^.il-'-^^m-"''. (87) 



X x^ 



Using this in (84) we get 



S + 2/2ir; = o, 

X ^ 



which on substituting for fi and /2 leads to 



X^P . fil^l^^ 



«(^A;.n(^ii^AVj^O. (89) 
Making use of the trigonometric half-angle formulae the above equation becomes 

where in writing the above equation we have ignored the possible sin((7^p)Pi^'^/2x^^) = 
solution since, as noted earlier, this corresponds to a trivial solution from the point of view 
of quantum dynamics. To solve this equation we make the ansatz = 2x'^p sin~^ (r) / {■yipY 
and solve for r to find r = c/x*^^^^^-*/^, where c is an arbitrary constant. This implies that 

2x'^P { c \ 



(74)P Vx(4p+i)/2 

To fix the constant c, we first note that it has dimension and secondly that 

when quantum effects tend to zero (as characterized by tending to zero for a fixed 
X ^ Vt^)? should tend to the classical value of zero. Then the obvious choice 
for c that suggests itself is c = {'jip)^^^'^^^^^. We thus have the complete solution with 
= —Nf2 after substituting for N, f2 and evaluating to 

V x^P+i /V a; x^ J ^ ' 

and we note that in the absence of quantum effects this correctly gives A^^ = 0. 

With the solution for the phase space variables in hand the obvious thing to do would 
be to write a metric corresponding to this solution. However, as seen in the last section for 
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version II of inverse triad corrections, with modified constraint algebra conventional space- 
time notions including the usual notion of metric do not make sense since the covariance 
property as encoded in the constraint algebra is modified. Nevertheless there certainly 
is gauge invariance in the theory since the constraint algebra is not anomalous, it is just 
deformed. Thus once again we need to identify what spacetime notions and covariance 
properties does the deformed algebra for the holonomy corrections imply. 

To aid in such an identification it would be convenient to have an alternative solution for 
the above set of constraints and the equations of motion since classically two such solutions 
are equivalent (being related to each other by general covariance). We therefore obtain 
another time independent solution which classically would correspond to the Painleve- 
Gullstrand gauge. For this we make the choice N = 1 which leads to E'^ = x and 



= —f2- Substituting for in (84) we get 



^ + 2/2if; = ^, (93) 

X ^ X'^ 



This equation is similar to (88) except that the r.h.s. in the present case is non-zero. We 



can solve this equation following the previous method to find 

- (tCJ)" 4x*+2 ) ■ * ' 

Appealing to the classical solution whereby K^p = —{2GM/x — GQ'^ /x'^Y^'^, we choose the 
minus sign in the above expression and fix ci = GM{'yip)'^P /2. 
given by 

f2GM GQ^^ ^''^ r 



\ X X 



2 



■ _ GM{ntlf^ GQ2(y|)2pn 1/2 



(95) 



which gives the correct classical limit = {2GM/x — GQ"^ / x'^y/'^ . 

With two solutions in hand the only task now is to find a suitable mapping from one 
to the other. However, for the present case this task appears to be much more non-trivial 
than for the inverse triad corrections. This can be understood by noting that presently we 



have only (84) to solve for K^p and A^. For any time independent choice for the form of say 
A^, we can solve for K^^ and it will be a valid time independent solution for the constraints 
and the equations of motion. This would also imply that for a given classical solution there 
are an infinite number of holonomy corrected solutions (in the sense that they all have the 
correct classical limit). 

For instance, above we found a solution for the choice A^ = 1 and noted that in the 
classical limit this would correspond to the Painleve-Gullstrand solution. But we can as 
well choose A^ = cos{{^^1,YK^p/x'^^) and to be specific we choose p = 1/2. In the classical 
limit a/ "yip I X — )■ and this choice will also go over to the Painleve-Gullstrand guage 
A^ = 1 (unless the solution of ( [84) ) for K^p has some singular behavior in that limit). For 
the choice A^ = cos( a/ '"yl^K^/x) however no such singular behavior occurs and the classical 
limit infact gives A^ = 1. 
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Thus we have the situation where corresponding to a given classical solution, there 
are (presumably) an infinite number of holonomy corrected solutions. It is therefore not 
clear that how or which two quantum solutions corresponding to two different categories 
of classical solutions (equivalently two different classical gauge choices) should be related 
to each other. To be specific, since both the choices = 1 and = cos{-\/'^lpK^p/ x) have 
the same classical behavior (and there would be many more such choices) it is not at all 
easy to decide as to which of these solutions should be related to the solution in equations 



(86)- (92) (with p = 1/2 for the case under discussion) which classically would correspond 
to the Reissner-Nordstrom solution in Schwarzschild-like coordinates. 

By the same token there would be many quantum corrected solutions which would 
correspond to the classical Reissner-Nordstrom solution in Schwarzschild-like coordinates. 
It is likely that with more quantum corrections included simultaneously, say the inverse 
triad corrections, this degeneracy will get broken and one will recover unique correspon- 
dence between the quantum and the classical solutions. Compared to their homogeneous 
counterparts where inverse triad corrections are found to be sub-dominant, this further 
highlights the significance of inverse triad corrections in the context of inhomogeneous 
systems though in a slightly different context. 

In the absence of suitable notion of a metric, it is unfortunately not feasible to explore 
the effecs of holonomy corrections on the properties of classical Reissner-Nordstrom black 
hole. We however note that as pointed out in [26l [3] , with holonomy corrections there is 



a possibility of signature change. Indeed, with /2 as given in (74), the deformation factor 



df2/dK^m (72) is 



For the solution in the Schwarzschild-like coordinates, for instance, using (91) this leads to 

dh , , 2(74)4p+i/2 



(97) 



which is negative for x < 2^^^P~^^{'yipY^'^ and would lead to Euclidean signature in the 
deep quantum regime. Similar conclusion will hold for other gauge choices also (signature 
change being a gauge independent concept), though in the absence of a suitable notion of 
covariance two different gauges cannot be directly compared. 



5 Conclusions 

We considered the effect of loop quantum gravity corrections on the properties of Reissner- 
Nordstrom black holes. Effects of both inverse triad corrections and point holonomy cor- 
rections were considered. It turned out that for both kinds of corrections an anomaly-free 
version of constraint algebra is possible. While for inverse triad corrections, version keeping 
the constraint algebra unmodified can be found, for holonomy corrections considered here 
this does not seem to be a possibility. As a further contrast between these two kinds of 
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corrections we found that for inverse triad corrected constraints and equations of motions, 
static solutions are possible generically (that is, independent of the specific form of correc- 
tion functions a and a) whereas for holonomy corrections static solutions do not exist and 
only stationary solutions are allowed. 

As has already been shown in [T5j, we once more confirmed (in the presence of Maxwell 
field) that for the version of inverse triad corrections resulting in a modified constraint 
algebra, spacetime covariance property is modified though this is not the case where only 
the dynamics is modified and not the algebra. It is a non-trivial task to identify what 
a modified version of spacetime covariance could be so that one could meaningfully talk 
about black holes and their horizons. For inverse triad corrections we were lucky and 
could identify such a notion based on quantum version of mapping from phase space to 
spacetime. Since in the quantum theory there is no spacetime to begin with, a possibility 
where the classical mapping emerges only in the classical limit cannot be ruled out. As 
required, this quantum version of the mapping goes over to the classical mapping when 
the quantum gravity effects are absent. 

For the case of holonomy corrected constraints, however, we could not find a modi- 
fied notion of spacetime covariance. A possible source for this difficulty was pointed out 
by showing that corresponding to a classical solution in a given gauge there correspond 
(possibly) infinite quantum corrected solutions which have the correct classical limit. It 
therefore becomes highly non-trivial to relate quantum solutions in two different gauges. 

However, it is possible that such a large degeneracy is present because we ignored 
other corrections. It is very much possible that including inverse triad corrections simul- 
taneously with holonomy corrections would lift the degeneracy. This further highlights 
the significance of inverse triad corrections in inhomogeneous situations compared to their 
cosmological counterparts. Unfortunately, for this reason we could not write a spacetime 
metric for the phase space solutions found for holonomy corrected constraints. However 
we found that as recently suggested in the literature, there exists possibihty of signature 
change with holonomy corrections. 

For both versions of inverse triad corrections, however, we could write a covariant metric 
and analyze some of its properties. We found that while for the case of unmodified algebra 
quantum gravity corrections resulted in mass, charge and wavefunction renormalization, 
for the case with modified constraint algebra only the charge was renormalized. A possible 
wavefunction renormalization (as suggested by phase space solution for the lapse function) 
did not appear in the metric because of the modified mapping used to go from phase space 
to spacetime. 

In both the cases it turned out that Newton's constant was not renormalized. Such 
a conclusion would have been difficult to draw from the quantum version of vacuum 
Schwarzschild solution since there the mass is multiplied with Newton's constant and a 
priori (at the effective level) it cannot be decided whether the correction function is mass 
renormalization or is it the renormalization of Newton's constant. Presence of Maxwell 
field breaks the conundrum and shows that G is not renormalized. We also found that 
inverse triad corrections lead to the possibility of three horizons over a finite range of mass 
in deep quantum regime and complete disappearance of the inner horizon (one of the dis- 
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tinguishing features of classical Reissner- Nordstrom blackholes) beyond a certain value of 
the mass parameter. 
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